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Abstract 

High temperature expansions for the susceptibility and the second correlation 
moment of the classical N- vector model (also known as the O(N) symmetric 
Heisenberg classical spin model) on the square lattice are extended from order 
/3 14 to (3 21 for arbitrary N. For the second field derivative of the susceptibility 
the series expansion is extended from order /3 14 to /3 17 . 

For —2 < N < 2, a numerical analysis of the series is performed in order 
to compare the critical exponents 7(iV), v(N) and A(iV) to exact (though 
nonrigorous) formulas and to compute the " dimensionless four point coupling 
constant" g r (N). For N > 2, we present a study of the analiticity properties 
of X) £ etc. in the complex f3 plane and describe a method to estimate the 
parameters which characterize their low-temperature behaviors. We compare 
our series estimates to the predictions of the perturbative renormalization 
group theory, to exact (but nonrigorous or conjectured) formulas and to the 
results of the 1/N expansion, always finding a good agreement. 

PACS numbers: 05.50+q, 64.60.Cn, 75.10.Hk 
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I. INTRODUCTION 



We have extended the high temperature (HT) series expansion of the zero field sus- 
ceptibility x{N',P) and of the second correlation moment /j, 2 {N;P) to order (3 21 and the 
expansion of the second field derivative of the susceptibility X4(iV; 0) to order f3 17 for the N- 
vector model Jl| (also known as the O(N) symmetric Heisenberg classical spin model) on all 
bipartite lattices in d = 2, 3, 4, 5, .. space dimensions The series coefficients have been 

determined by using the vertex renormalized linked cluster expansion (LCE) method 0-0] 
and have been expressed as explicit functions of the spin dimensionality N. This calculation 
pursues and improves our previous work to a considerable extent: it summarizes into a 
convenient format a large body of information for an infinite set of universality classes and 
offers further insights into the properties of the iV-vector model by enabling us to vary with 
continuity the crucial parameter N and to study how various physical quantities depend on 
N. 

Strictly speaking, the N— vector model is defined only for positive integer N. Therefore 
it is possible to construct infinitely many "analytic interpolations" in the variable N of the 
HT coefficients and, as a consequence, of the physical quantities. We have performed the 
"natural" analytic interpolation by which the HT coefficients are rational functions of N. 
This is the most interesting interpolation because it coincides with that used in the 1/N 
expansion as well as in the usual renormalization group (RG) treatments and moreover it is 
unique in the sense of Carlson theorem |Ttj . 

Next interesting step, on which we will report elsewhere [|J, is to compile tables of HT 
coefficients analytically interpolated both in N and in the space dimension d. The "natural" 
analytic interpolation, with respect to d, of the HT coefficients, which is polynomial in d 
and equivalent [|ll|] to the one of the Fisher- Wilson method []12 |, is also unique in the sense 
above specified. We will thus be able to describe accurately the general (N, d) universality 
class, to achieve also a "view from HT" of presently inevitable RG approximation schemes 
such as the e expansions (at the upper and at the lower critical dimension), the 1/N and 
the 1/d expansions and possibly to gain a more detailed knowledge of their limitations. As 
pointed out in Ref. [0,0 possible violations of the convexity of the free energy or of the 
Lee- Yang property might occur for noninteger values of N or for noninteger values of d, 
respectively, but, of course, they do not question the long known conceptual and analytical 
advantages of treating iV or d as continuous parameters. 

This paper is devoted to the study of the iV-vector model on the square lattice and 
its main result reported in the appendices, is the tabulation of the HT series expansion 
coefficients. We have made an effort to keep our exposition selfcontained, in particular 
by reporting all known HT expansion coefficients of \i A*2 and %4 an d not only the newly 
computed ones. It is worth noting that further sizable extensions of these series are not 
too difficult and now are in progress. Several interesting, but somewhat intricate, computa- 
tional procedures that have made this laborious calculation (and its forthcoming extensions) 
possible will be illustrated elsewhere ||. 

We shall also update, refine and extend our previous numerical analysis of 0((3 1A ) series 
presented in Ref. ||, but our discussion will be more sketchy whenever there is some overlap 
with that reference. 

The paper is organized as follows. In section II we introduce our notations and definitions. 
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In section III and IV we present a numerical study of our new HT series which are now long 
enough for a reliable assessment of the uncertainties of the analysis initiated in ||||. 

In section III we examine the range —2 < N < 2 in which the iV-vector model is known 
(or expected) to have an ordinary power law critical point at finite temperature and we 
estimate the critical exponents j(N) of the susceptibility, v(N) of the correlation length, as 
well as the gap exponent A(iV) from Xi i n order to compare them with exact (nonrigorous) 
formulas proposed some time ago ]I3]-|nj. We do not reanalyse here the N = 2 case (the 
Kosterlitz-Thouless model) which has been already studied in Ref. [jD]], but we only wish to 



call the readers attention on a few tiny numerical errors which have crept in the HT series 



coefficients at orders 17-20 reported in [I7J and were due to an accidental contamination of 
a numerical file in the final stage of that calculation. In this paper we have corrected such 
errors, which of course could be annoying to those who wish to extend the computation and 
first have to make sure that they are able to reproduce correctly the existing data. We have 
also checked that, being so small, these errors were of no consequence at all either on the 
qualitative and on the quantitative results of the analysis in Refs. jTTj-fZIJ, which therefore 
does not need to be repeated until significantly longer series and/or better methods of 



analysis can really offer new insights, for example on the questions raised by Ref. |22] 



In section IV we examine the set of models with N > 2, which are expected to behave 
quite differently ||23|| . In the last two decades their features have been extensively explored by 
various analytical and numerical techniques, with the main motivation that they are lattice 
regularizations of the field theoretic non-linear 0(N)-symmetric cx-models, which share the 
crucial asymptotic freedom property with four dimensional gauge field theories, but are 
much easier to study. We update our previous survey of the nearby singularities of Xi ^2 
(and Xi) m the complex inverse temperature (3 plane || and we still find no indication 
of any physical critical point at finite real (3. On the contrary, we point out that the low 
temperature behavior appears to join smoothly onto the high temperature behavior so that 
several parameters which characterize the low temperature behavior can be computed in 
terms of HT series and full consistency is obtained with the predictions of the perturbative 
RG. 

We end the paper by comparing our conclusions to some related recent works which, 
either by direct stochastic simulations or by analytic approximations such as the 1/N ex- 
pansion, also test and confirm the predictions of RG. 

In the appendices we report the closed form expressions for the HT series coefficients of 
X, H2 and Xi as functions of the spin dimensionality N and their evaluation for a few specific 
values of N. Electronic files containing these data may be requested from the authors. The 
present tabulation extends and supersedes the one in Ref. M which, unfortunately, is marred 
by a few misprints. 

II. DEFINITIONS AND NOTATIONS 

The Hamiltonian H of the iV-vector model is: 

H{v} = - 1 - £ v(x).v(x'). (1) 

(x,x') 
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where v(x) is a iV-component classical spin of unit length at the lattice site x, and the 
sum extends to all nearest neighbor pairs of sites. 

The HT expansion coefficients of all correlation functions C(N;/3) = J2 r fr{N)f3 r are 
simple rational functions of N, namely f r (N) = P r (N)/Q r (N) where P r {N) and Q r {N) are 
integer coefficient polynomials in the variable N. Therefore the same property is true for 
the expansion coefficients of the susceptibility 

oo 

X {N, (3) = 5>(0) • v{x)) c = 1 + E ar(N)[3 r ; (2) 

x r =l 

of the second correlation moment 

oo 

» 2 (N,(3) = $»(0) • v(x)) c = J2 s r (N)p r (3) 

x r=l 

of the second field derivative of the susceptibility 

1 oo 

X*(N, (3) = £ (v(0) ■ v(x)v(y) ■ v(z)) c = - ( - 2 + £ d r (N)(3 r ) . (4) 

x,y,z iv r =l 

It should be noticed that our definitions of x and fi 2 differ by a factor and our definition 
of Xi differs by a factor N ^ +2 ^ from that of Ref. (apart from misprints, also adopted in 
Ref. || as far as %4 is concerned). 

Like in any calculation of HT series, the correctness of the numerical results is a decisive 
issue. Our confidence on the validity of this work is based not only on the numerous direct 
and indirect internal tests passed by our code, but especially on the fact that the space 
dimension d and the spin dimensionality N enter simply as parameters which can be varied 
to produce, always by the same procedure, series in complete agreement with those already 
available (sometimes to a higher order) for specific values of iV and d and in the spherical 
model limit (namely the limit as iV — > oo at fixed (3 = (3/N) More precisely let 
us observe that from the tabulation in the Appendix it appears that d r (N) = N r a r (N), 
s r (N) = N r s r (N), and d r (N) = N r d r (N) have a finite limit as iV — > oo, so that the spherical 
model limit of our series coefficients can be immediately read. The HT coefficients of the 
susceptibility x^{0) an d the second correlation moment ji^ {f3) of the spherical model can 
then be used to check in our tables the coefficients of the highest power of iV in the numerator 
polynomials of a r (N) and s r (N) respectively. It should be noticed that %4 is 0(1/N) for 
iV — >• oo, as it is expected because, in the spherical model limit, only the two-spin connected 
correlation functions are nonvanishing. However the quantity xa(N,/3) = Nxi{N,f3) has a 
finite large N limit which will be denoted by x£ (fl)- We remind the interested reader that, 
in Ref. [23], we have tabulated the HT coefficients for the spherical model susceptibility x^ 
on the square lattice through order 63, and that the HT coefficients for yffi and xi can 
be obtained from the expansion of x by using the formulas n { 2 s \{3) = A/3( X {s) 0)) 2 and 

xi s) 0) = -2( X ^0)) 2 (x {s) 0)+^))- 

Similarly also the N = limit, which corresponds to the self-avoiding walk (SAW) model 
can be easily obtained from our HT series after expressing all quantities in terms of $. Again, 
it is the quantity X4(N, (3) that has a finite limit for iV = 0. The HT coefficients are then 
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essentially given by the constant terms in the numerators of a r (N), s r (N), and d r (N), up 
to simple factors from the denominators. 

More (partial) checks are trivially obtained by setting AT = 1, 2 or 3 in our formulas and 
comparing the results to the corresponding available expansions cited below. 

Of course, if a "complete set" of such checks were available, it could be used to reconstruct 
the whole computation. 

It is interesting to recall that, more than two decades ago, HT series valid for all iV and 
for a general lattice have been computed up to order /3 8 ( j3 9 for loosely packed lattices) 
2l| . Later on, in the case of a square lattice, the series were extended through (3 11 || and 



then through /3 14 0,0]. On the other hand, very long expansions, on the square lattice, have 
been computed recently, both for the susceptibility and for the second correlation moment 
in the special cases iV = f27|j2§ (the self-avoiding walk (SAW) model), and N = 2 ( 



the Kosterlitz-Thouless model), by highly efficient algorithms, whose performance, however, 
does not excel for space dimensionality larger than two or which cannot be extended to other 
values of N. More precisely, for the susceptibility, the published series extend through orders 
/? 43 (recently pushed to /3 51 |[29|| ) for A^ = and ft 20 for N = 2, and for the second correlation 
moment through /3 27 and j3 20 , respectively. The longest published expansions of valid 
for any N, presently extend only through j3 14 ||. The N = 1 case (the spin 1/2 Ising 
model), which is much simpler because the model is partially solved, should be considered 
separately: in this case the available series for x and \ii extending to /5 54 are tabulated in Ref. 
|30| while the series for X4 extends only to (3 17 . When our work was being completed, 



another calculation valid for any A" was announced |32[ for the nearest neighbor correlation 
function, x an d ^2 (but not for % 4 ) by the technique of group character expansion. This 
procedure seems to be efficient only in 2 space dimensions and to be presently feasible up 
to order 21 on the square lattice, to order 30 on the exagonal lattice and to order 15 on 
the triangular lattice f33|. It is reassuring that our general results, although obtained by a 
completely different procedure, agree throughout their common extent also with the specific 
cases A^ = 2, 3, 4, 8 tabulated in Ref. p3|. 



III. ANALYSIS OF THE HT SERIES FOR N < 2 

We will now discuss some of the information that can be extracted from the series and 
update the analysis first presented in |§. 

Let us recall that an exact expression for the critical exponent v(N) within the range 



—2 < A^ < 2 has been conjectured in Ref. [14] on the basis of an approximate RG analysis. 



Later on the same expression and an analogous one for the exponent r](N) were derived [T5| 
by observing that a special O(N) spin model (assumed to be a faithful representative of this 
universality class) can be mapped into a soluble []TB[ loop gas model. 
The conjectured exact exponents are 

HN) = ji^ (5) 

and 

j(N) = (2 - v (N))v(N) = (6) 
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with N = —2 cos(^) and 1 < t < 2. 

The quantities Xi £> an d X4 are then expected to display, in the whole range —2 < iV < 2, 
as (3 | Pc(N), the conventional power law critical behaviors x — c x (iV)(/3 c (iV)— /i) -7 ^, £ ~ 
Pe(iV)(/3 c (iV) - P)-^ N \ and x 4 ^ -c 4 (N)(f3 c (N) - p)-^)-2A(N)^ with Wegner » con fluent" 
corrections. Here c x (N), c^(N) and Czt(iV) are (nonuniversal) critical amplitudes. 

In order to test numerically the validity of @ and (|]), we have estimated 'y(N) and v(N) 
by forming first order inhomogeneous differential approximants (DA) of the susceptibility 
X and of the "second moment" correlation length squared £ 2 = H2/2dx respectively. We 
also have computed the gap exponent A(iV) from % 4 . Our numerical procedure consists 
in averaging over all estimates from DA's in the class selected by the protocol of analysis 



of Ref. |35[ which use at least 16 series coefficients in the case of x an d £, and at least 14 
coefficients in the case of X4- F° r each value of N, we first estimate (3 C {N) and 7(A) from x, 
then we use (3 C (N) to bias the computation of v(N) from £ 2 and of A(iV) from Xa/x- As a 
measure of the uncertainties we have taken three times the rms deviation of the approximant 
estimates. 

Alternatively, we have assumed the validity of the hyperscaling relation 

2A(A) = 2u(N) + j(N) (7) 

for —2 < N < 2 and used also the series for Xi 111 the computation of v(N) by re- 
sorting to the so called "critical point renormalization" In this case we have esti- 



mated v(N) by examining the singularity at z = 1 of the series J2r h r (N)z r with coefficients 
h r (N) = d r (N)/t r (N) where x 2 {N] (3) = E r U(N)p r . Similarly the exponent A(N) has been 
determined from the series with coefficients l r {N) = d r (N)/a r (N) and the exponent 'y(N) 
has been obtained in terms of the series for /12 and for £ 2 . This procedure does not require 
the knowledge of /3 C (N), but only seventeen term series are available for the computation of 
v(N) and A(N). 

In Fig.l we have reported our results for 'y(N), v{N) and 2A(iV) versus and compared 
them to the exact formulas (|^), (H) and to (0). 

In the central region of the plot, approximately for — 1 < N < 1.5, both numerical 
procedures we have followed yield very accurate estimates agreeing with the exact formulas 
within a small fraction of a percent. Near both ends of the interval —2 < N < 2 the 
agreement deteriorates because the series have to crossover to different singularity structures 
in order to exhibit either a gaussian behavior for iV = —2 or a Kosterlitz-Thouless behavior 
for N = 2. However, since the exponent estimates always move in the right direction as the 
number of series coefficients is increased, we are confident that, in these border regions, we 
are simply facing a numerical approximation problem rather than a breakdown of the exact 
formulas (|5|), (|J) and therefore we can conclude that their validity as well as the validity of 
the hyperscaling relation (0) are convincingly supported also by our HT series study. 

In the N = case our expansion for X4 is the longest presently available and therefore 
it is worthwhile to update the verification of the hyperscaling relation (|7|). If we bias the 
first order DA's of x^/x using the value /3 C (0) = 0.3790525(6), obtained in Ref. pSJ from 



an 0(/3 43 ) series for x, we get the estimate A(0) = 1.422(1) which is within 0.1% of the 
predicted value A(0) = 91/64 = 1.421875. Similarly (and with the same bias), from a 
study of Xilx 1 ! j we obtain the estimate: 2A(0) — 7(0) = 1.503(9) which by (0) and (g) 
should be compared to the exact value 2z/(0) = 1.5. By studying directly X4> we obtain 
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the estimate 2A(0) + 7(0) = 4.175(3). Adding the last two estimates, we can conclude that 
A(0) = 1.419(3), which is slightly less accurate, but perfectly compatible with the previous 
result. 

We can also estimate with fair accuracy the (nonuniversal) critical amplitudes of Xa 
for N = and N = 1 which might be useful for reference and comparison with other 
numerical calculations. Let us recall that in the Ising model case the critical amplitude of 
the susceptibility c x (l) has been computed exactly to be c x (l) = 0.962581732.. and that 
also the amplitudes of the first few subleading confluent corrections to scaling are known 
Since the first confluent corrections are found to be analytic, and indeed it has been 



argued |38] that there are no irrelevant-variable corrections to scaling in the thermodynamic 
quantities of the two-dimensional Ising model, we expect that we can rely quite simply on 
near diagonal Pade approximants (PA's) of {(3 C — /3) 7 ^x to obtain an accurate estimate 
of c x (l), even from not particularly long series. Indeed, from our 0((3 21 ) expansion, we 
get c x (l) = 0.96261(3) using as a bias the known values of 7(1) and /3 C (1). The critical 
amplitude of Xi has not yet been evaluated exactly, but since the structure of the confluent 
corrections to scaling should be similar to that of x also this amplitude should be accurately 
estimated. Our biased estimate of this quantity is £4(1) = 4.378(2) which compares well 
with the estimate £4(1) = 4.37(1) from the fifteen term series of Ref. ||39|| . Also in the N = 
case ( now in terms of the variable f3) the structure of the confluent corrections is likely 
to be favorable since both the long expansion computed for x i n Ref- and the results 



of extensive stochastic simulations |4(| are consistent with an analytic dominant confluent 
correction. In this case we get c x (0) = 1.0524(8) and £4(0) = 6.62(2). 

In terms of Xi £5 and Xi we can also compute the " dimensionless renormalized four point 
coupling constant" g r (N) as the value of 

at the critical point /3 C (N). If we assume that the inequality 7(iV) + 2u(N) — 2A(N) > 0, 
(rigorously proved to hold as an equality for N = 1), is also true for any —2 < N < 2 , 
then g r {N) is a bounded (nonnegative) universal amplitude combination whose vanishing is 
a sufficient condition for gaussian behavior at criticality, or, in lattice field theory language, 



for "triviality" of the continuum limit theory defined by the critical N— vector model [25 



Notice that our normalization of g r {N) is the same as the one adopted in Ref. and 
differs by a factor g7rJ ^J v 8 +2 ^ from the normalization traditionally chosen in the field theoretic 



renormalization group treatments |43| . 

For < iV < 2 we have estimated g r (N) by evaluating both near diagonal PA's and first 
order inhomogeneous DA's of the series for l/g r (N, (3) at the critical values (3 C {N). The two 
procedures yield results which are perfectly consistent within their numerical uncertainties. 
In Fig. 2 we have reported our estimates of g r {N) for various values of < iV < 00 and 
compared our results to other computations in the literature. 

For N = we estimate g r {0) = 10.53(2), which agrees well with a previous estimate 
&.(()) = 10.51(5) from the 0((3 U ) series 0|| studied in Ref. |4 . 



For N = 1, our estimate g r (l) = 14.693(4) coincides with the estimate of Ref. |5T 
obtained using the same number of coefficients, but a rather different method of analysis, and 



is consistent with previous estimates from an eleven term series |39[] giving ^ r (l) = 14.67(5) 
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and from a fourteen term series |44| giving g r (l) = 14.63(7). Less precise, but consistent 



estimates also come from field theory (fixed dimension RG) g5[ yielding g r (l) = 15.5(8) and 
from a recent single cluster MonteCarlo (MC) simulation El yielding g r (l) = 14(2).. 

For N = 2, biasing the approximant with the critical inverse temperature /3 C (2) = 
1.118(4), we get the estimate g r (2) = 18.3(2). This result is consistent both with the deter- 
mination g r {2) = 18.2(2) obtained in Ref. Q and with the MC measure g r (2) = 17.7(5) 
obtained in Ref. [0]. 

It should be noticed that, as a reflection of the growing complexity of the critical singu- 
larity structure, the uncertainty of our numerical results is very low for N = 1 and relatively 
modest for N = 0, but it is much larger in the N = 2 case. However our new estimates 
appear to be generally more accurate than previous ones. The corresponding calculation for 
N > 2 will be discussed in next section. 



IV. ANALYSIS OF THE HT SERIES FOR N > 2 



In this range of values of N the general features of iV-vector model are expected to 

D, 



change qualitatively: reliable, although nonrigorous (and sometimes questioned p3| . 
RG calculations at low temperature 
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indicate that the model is asymptotically free, 
namely that it becomes critical only at zero temperature. The asymptotic behaviors of the 
" second moment" correlation length, of the susceptibility and of xa as P — ► oo are predicted 
to be 



r(N,P)=c^N)( 



-p . 

bo(N)- 



exp\ 



~P - 

b (Ny 



i + 



Hi(N) 

P 



p 2 K p*' 



(9) 



X as {N,(3) = c x {N){ 



-P ■ 
b (N)' 



a»i(JV) 

) b (N)2 



b (N) 



exp\ 



-2P_ 
b (N) 1 



Kt(N) K 2 (N) 



P 



P 2 



(10) 



xt(N,P) = -c 4 (N)( 



-3 6blW I 2g (JV) 



bo(N) 



exp\ 



where 



bo(N) 



jN-2) 
2tt 



h{N) 



<N-2) 
(2vr) 2 * 



-6/3 
6 (iV) J 



9o(N) 



P' 



N 



2tt 



(11) 



(12) 



are the first (renormalization scheme independent) low temperature perturbation expan- 
sion coefficients of the RG beta and gamma functions [|23[ and c^(N), c x (N) and c 4 (iV) are 
universal quantities which clearly cannot be computed in (low temperature) perturbation 
theory. The (nonuniversal) constants Hi, H 2 , Ki and K 2 appearing in ([]) and QTD|) can be 
calculated in low temperature perturbation theory, and, on the square lattice EJ^j, they 
come out rather small but not completely negligible in the range of values of /3 in which we 
shall be able to compute reliably x an d £ • Numerical estimates for Hi, H 2 , Ki and K 2 can 
be found in Ref. and for brevity are not reported here, although we use them in the 
calculations. Unfortunately, the analogous 0(4) corrections have not yet been computed for 
Xa- As a consequence of (0) and flTTD , for large P, 
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m) = ^ (J v)[i + 0( i)] = ^g^ n + o ( i)]. us) 

where c 4 (iV) = Nc A (N). 

Let us notice that the asymptotic formula (|9|) is valid both for the "second moment" cor- 
relation length £(N, (3) (which here is employed exclusively), and for "exponential correlation 
length" £ exp (N,j3) but, a priori, with different multiplicative constants c%(N) and c e ^ xp (N), 
respectively. However, it has been repeatedly noticed that £ exp and £ are numerically very 
close in the critical region ||50|| . This fact is confirmed by a recent analytic calculation for 
large iV of the (universal) ratio c e ^ xp '(N) / 'c^(N) giving the result | 5T[ 



cf p (N)/c^N) = 1 + 0.003225/iV + 0(1/ N 2 ) (14) 



Moreover, for N = 3, this ratio has been measured [0] by a high precision MC method 



at /3 — 1.7 and 1.8, fully confirming the quantitative reliability of the 1/N expansion flU]) 
down to very low values of N . Therefore, with very good approximation, we are justified 
in simply identifying c^ xp (N) and cg(iV) even for small N. 

These results are of direct interest here because the coefficient c^ xp (N) can be computed 
exactly by the thermodynamical Bethe Ansatz and its value, with our normalization 
conventions, is 



r(i + 



c 



T P (N) = ^- 2 ' exp - '-f^ll ~ c^ N ) (15) 



32 



n/2 + log(8/e)) 



N-2 



Let us now turn to series analysis. With sufficiently long series, like those analyzed 
here, even the simple plot of the HT coefficients of Xi versus their order, reported in Fig.3, 
is suggestive enough that the series for iV > 2 must behave quite differently from those 
for N < 2. Notice that in the plot the coefficients have been conveniently normalized, 
as indicated in the figure caption, in order to make the behaviors for different values of 
N easily comparable. In the N < 2 case the coefficients are positive (and monotonically 
increasing) in agreement with the fact that the nearest singularity is located on the real 
positive (3 axis (and that the antiferromagnetic singularity located at —/3 C is much weaker). 
On the contrary, in the N > 2 case the coefficients do not remain positive as their order 
grows and they display what can be safely interpreted as the onset of an oscillatory trend. 
This feature is related via Darboux theorem to the fact that, for N > 2, the nearest 
singularities of x i n t ne complex (3 plane become unphysical, as we have first pointed out 
some time ago 0. This applies also to /i 2 , as well as to Xi with the only difference that for 
these quantities the oscillating behavior of the expansion coefficients should set in at higher 



orders. Our interpretation of these general features is impressively confirmed by a study [24 
of the spherical model which, in spite of a priori legitimate mistrust about exchanging the 
large (3 and the large N limits, turns out to be a completely reliable guide to the qualitative 
behavior of the iV— vector model even for not too large iV > 3, as it has been also argued 
some time ago [p|,|55 



As already mentioned above, an arbitrarily large number of HT expansion coefficients 
for x^ s \ a4 an d Xi °f the spherical model can be easily computed, for any lattice. In the 
case of the square lattice they exhibit regular cyclic alternations in sign of period 8 related 
to the symmetric quartet structure of the nearest unphysical singularities in the complex 
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(3 plane. We have accurately mapped out in Ref. p4j the whole set of singularities, all of 
which are square root branch points. In the vicinity of (3 = oo, this set has the characteristic 
structure dictated by asymptotic freedom (which was first discussed for the case of QCD 
in Ref. J56|), the analiticity domain of x > being a wedge with zero opening angle which 
contains the real j3 axis. It is quite likely that these features of the spherical model persist 
also down to all finite N > 3 ]57|], although a complete study of this question is presently 
infeasible. 

The transition from the N < 2 regime characterized by a conventional power law critical 
point, to the N > 2 asymptotically free regime characterized by unphysical singularities can 
be closely followed by locating the position of the nearest singularity (3 n in the first quadrant 
of the complex (3 plane as a function of N. For convenience in the graphical representation of 
the results, we shall use in what follows the scaled variable (3 = (3/N and plot the estimates 
of (3 n versus x = 1 — 1/N, rather than versus N. The trajectory of the singular point f3 n as a 
function of x in the complex (3 plane can be traced out as described in Ref. || either by using 
PAs to locate the nearest singularity of the log derivative of x or by directly computing DA's 
of x- The results of both procedures agree perfectly within the numerical uncertainties. 

In Fig. 4 we have plotted the real and the imaginary part of (3 n as functions of x in the 
range < x < 1. For —2 < N < N < 2.2, the singularity (3 n is still a real critical point, 
but for iV > N it splits into a pair of complex conjugate singularities which move into the 
complex plane and, as iV — ► oo, tend to the limiting points (3± ~ 0.32162(l±z). In particular 
for iV = 3 the nearest pair is located at (3 n = 0.58(5) ± 0.14(5)z, while for N = 4 we have 
(3 n = 0.55(4) ± 0.22(5)z, for iV = 6 we have (3 n = 0.50(2) ± 0.28(2)i and for iV = 10 we 
have j3 n = 0.44(2) ± 0.31(l)z. Although in general they may be weaker, the corresponding 
antiferromagnetic singularities will follow trajectories symmetrical with respect to the Im(f3) 
axis so that, for all iV > N the set of the nearest singularities will form a quartet with the 
same symmetry. 

It is certainly conceivable [4~B] that, in contrast with the perturbative RG predictions, 
when the nearest singularities become complex, a further real critical singularity might 
appear so that, even for iV > 3, it would be still possible to relate the steep dependence 
on (3 of Xi £ an d Xi t° a conventional finite temperature phase transition, but we have not 
been able to find any numerical indication of such a possibility for not too large j3. More 
precisely neither Dlog PA's nor DA's exhibit any real and numerically stable singularity 
in their range of sensitivity. Another argument against the existence of critical points for 
finite values of (3 comes from the observation illustrated below, that, by a simple procedure, 
the high temperature behavior of £, x an d Xa can be smoothly extrapolated onto their low 
temperature behavior (0), (PD and flTTD as predicted by the perturbative RG. This is feasible 
for any N, although the procedure is numerically very accurate only for N > 3, since, for 
N = 3, Im((3 n ) is small and therefore the behavior of x or £ on the real (3 axis is more 
strongly perturbed in the vicinity of Re((3 n ). Both the often reported failure in observing 
asymptotic scaling by MC simulations of the N = 3 model at moderate values of (3 and 
the better successes for larger values of N find a completely plausible explanation in this 
picture. Of course, the results of our extrapolation scheme would be difficult to explain if 
the high temperature region were separated by a critical point from the low temperature 
region. 

Let us now describe an approximation scheme which enables us to estimate low tem- 
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perature perturbative parameters such as bo(N), as well as nonperturbative parameters like 
c^(N), c x (N) etc. entering into the asymptotic formulas @, (|T0| ) and ([TTJ) , in terms of our 



HT series. Since £, x etc. are exponentially fast varying quantities at large values of j3, 
neither PA's nor DA's are well suited for a straightforward extrapolation of the HT series 
from small to (relatively) large (3 values. We should rather work with quantities which vary 
slowly enough to be well represented by PA's or DA's. Let us observe that, if x has the 
asymptotic behavior fllPp, then for large enough /3 

B X (N, (3) = l -Dln[ x {N, 0)] + - ±Dln[l + ^| + -§-] * -N/b (N) + 0(- ' 



2(N-2)P 2 N(3 N(3 2 iV 4 /? 4 

The log derivative of Xi which is a slowly varying quantity, can be approximated by near 
diagonal PA's and then we can reliably extrapolate the quantity B X (N,(3). In practice, due 
to the finite extension of our series and to the intricate analytic structure of Xi we do not 
expect that this is a good approximation for large (3 and we rather make the reasonable 
(and successful) assumption that the f3 independent parameter b (N) is best approximated 
by evaluating B x (N,/3) at some finite real value (3 = (3 S where it is stationary or it shows 
the slowest variation when (3 is varied. Consistency of this approximation scheme requires 
that, as the number of HT coefficients used in the calculation is increased, the stationary 
value B X (N,/3 S ) stabilizes and that (3 S — > oo. It can be checked that this actually happens 
in the N = oo case in which arbitrarily long HT expansions can be studied and also that 
our approximation scheme converges rapidly to the expected result. 

A further check of the correctness of our procedures comes from the obvious remark that 
similar estimates of the same parameter b Q (N) should be obtained starting either with the 
correlation length £ and computing the quantity 

B^N, ft = l-Dlnim pflP] + — J* - Dln[l + -5, + -f§-] ~ -N/b (N) + 0{ ' 



2 2(N-2)/3 N(3 N 2 f3 2 N 4 [3 4 

or starting with % 4 and computing the quantity 

Bi(N,~P) = \Dln[ X ,{NM + ^ ~N/b (N) + 0(j^). (18) 

In Fig. 5 we have plotted B X (N,(3) versus (3 for various values of N, in order to show 
that a stationary point f3 s (N) actually exists around (3 ~ 0.55 for any N, and that the size 
of the neighbourhood of $ S (N) where B^(N,(3) varies slowly with /3, grows with N. Notice 
that Re(f3 s (N)) < 0.55 for > 4 and therefore on the border of the convergence region of 
the series or slightly outside it. 

In Fig. 6 we have plotted versus iV our numerical estimates of bo(N) from the quantities 
B x (N,j3 a (N)),B((N,j3 a (N)) and B 4 (N,(3 S (N)) and have compared them to the expected 
value (0). Each point represents the average of the near diagonal PA's using at least 14 
series coefficients for the quantities B x (N,f3 s ) and B^(N,(3 S ), and at least 10 coefficients for 
B^{N, (3 S ). We have reported relative errors of 5% suggestive both of the scatter of the esti- 
mates obtained by the various PA's and of the systematic uncertainties of our extrapolation 
procedure. 

In conclusion, it appears that from our high temperature data for £, x an d X4? we have 
been able to extract completely consistent and correct estimates of the low temperature 
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perturbation parameter bo(N) which characterizes the exponential asymptotic growth of 
these quantities, and in general that the deviation from the expected value (0) of b (N) is 
never larger than 5% over a wide range of values of N. 

In quite a similar way, assuming that bo(N) is given by (|12"D , we can estimate the expo- 
nents of the power law prefactors in (^), (|TCi| ) and ([[J]). As it must be expected, the errors 
in this computation are somewhat higher, but they do not exceed 20 — 30%. 

Let us now show that by a similar approximation procedure we can also estimate the 
constant c^(N). We have simply to compute the HT series of the slowly varying quantity 



C,(N,P) = e(Nj)exp[^}^ (19) 
, ^,-N(3,^ f H x H 2 1 



obtained by dividing out the exponential factor in the asymptotic behavior (|jD of £ 2 . 
We then form near diagonal PA's to C^(N, $) and use them to evaluate the quantity 

C,(Nj)(^l)^(l + + J^Y 2 (20) 

at the value 3 S where it is stationary. In analogy with the previous computation this 
is taken to be the best approximation of cg(iV) 2 . It is observed that also in this case the 
stationary values occur for (3 ~ 0.5. 

Similarly, we can estimate c x (N) by studying the HT series for the quantity 



C x (Nj)= X (Nj)exp[^-}^ (21) 



-N3 -N-i / Ki K 9 1 



Unfortunately no exact formula is known for c x (N) , but we can compare our numerical 
estimates to the 1/N expansion through 0(1/N) of c x (N) which has been computed in Ref. 

m 



Cx(A0 = ^[l-^ + O(l/iV 2 )] (22) 



or to an analytic formula recently guessed in Ref. ||BT]] with no other theoretical justifi- 
cation than a formal analogy with the exact formula (|T5| ) for c^(N). 
An estimate of c 4 (iV) could be obtained starting with the series for 

6N3 —N3 -2(n+2) i 

c l{ N, n e^m^i^jo^)— (23) 

However the O(jj^) corrections are not known, and moreover the 1/N expansion of c^N) 



which has been computed in Ref. E| 
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(24) 



it is practically useless, except for very large N, since the subleading term is quite large. 
Therefore we do not report our estimates for c^N). 

By the same method we have also directly estimated the universal quantity c r (N) = 
^c^(N)/c x (N) which appears in the asymptotic expression of the ratio 



-N-2 

-2 



L * c r (N)^(l + ^ + -i + 0(^)) 2 (1 + *± + -*L.)~\ (25) 

X Nf3 N 2 f3 2 ( N 3 f3 3 N{3 N 2 f3 2 

and have compared it to the 1/N expansion 

c r (N) = 1 + 1.955/iV + 0(1/N 2 ). (26) 

Let us notice that, for large N, the unphysical singularities of £ 2 and x tend to cancel in 
the ratio and that the 1/N correction in ( p6|) is not very large. 

We have reported our numerical estimates for c%(N), c x (N) and c r (N) in Table 1 and 
in Fig. 7 where they are compared to the exact or conjectured formulas and to their 1/N 
expansions. Like in the previous Fig.6, the error bars we have attached to our data points 
are fairly subjective in that they include a "statistical" contribution (describing the spread 
of the estimates from various approximants) which is not large in general, while the main 
part of the uncertainty comes from our estimate of the systematic error. As it appears from 
the Fig. 7 and from Table 1, our central estimates for c%(N) and the exact formula agree 
within 1 — 2% on the whole range of N except for the lowest values of N. It should be 
observed that we have not our reported estimates for N = 3 since in this case the nearby 
unphysical singularities have a very small imaginary part Im(/3 n ) and there is a large spread 
in the stationary values of C% and C x . This makes difficult to estimate unambiguously the 
values of and c x . However, if we shift to only slightly larger values of N, such as iV = 3.5, 
then Im(p n ) is already sufficiently large for our procedure to work appropriately and we 
can estimate cg(3.5) = 0.028(8) to be compared to the exact value c| xp (3.5) = 0.0273 and, 
similarly, c x (3.5) = 0.021(2), while the conjectured formula gives c x xact (3.5) = 0.0266. In 
both cases the discrepancy is less than 20%. At iV = 4 the exact value is c| xac *(4) = 0.0416 
and we find c^(4) = 0.039(1), which is off only by 6%. For larger N the agreement is much 
closer as it is shown in Table 1. Our estimates also agree well with the conjectured exact 
formulas for c x (N) and c r (N). The discrepancy from these formulas or from their 1/N 
expansions does not exceed 5 — 10% for small values of iV but it gets significantly smaller 
already for moderately large N. All numerical results are collected in Table 1. We believe 
that both the failure to reproduce accurately the N = 3 values of the parameters and the 
other general features of our approximations should not be surprising if we take into account 
the analytic structure in the f3 complex plane of the quantities to be extrapolated and we 
consider that our computational method is the simplest and most direct possible and also 
that we are still working at moderate values of (3 where, for small values of N, the correlation 
length is not very large. Our approximation procedures should not however be suspected to 
be " ad hoc" , since they were proposed and the first results |9| were published before either 
the exact formula (O) and the 1/N expansions became known. 
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We should also at this point recall our remark |24| (resumed and applied to fourteen 
term series in ||62|| ) that the precision of these estimates might be significantly improved 
by performing a conformal transformation of the complex {3 plane in order to remove at 
least the quartet of the nearest unphysical singularities before applying our approximation 
procedure. We think, however, that the success of our straightforward treatment over a wide 
range of values of N cannot be accidental and that simply getting a higher level of accuracy 
could hardly be more convincing of the validity of the RG picture of scaling and indicative 
of the purely numerical origin of the discrepancies for the lowest values of N. Therefore we 
shall not pursue here our old suggestion. 

Finally, we can also estimate g r (N) for N > 3, by forming PAs to the series expansion 
of l/g r (N, f3) and evaluating them at their stationary points. We have reported in Fig. 2 our 
estimates and have included for comparison the field theoretic estimate for N = 3, and 
MC estimates |63|j63j42l| and other HT series estimates E3 for N = 3 and iV = 4. It should 



also be noticed that our results are entirely consistent with the large N limit, in which we 
have 



g r (N) = 8tt[1 - 0.602033/iV + 0(1/N 2 )]. 



(27) 

The 1/N correction has been computed recently ||44|| . Also the accuracy of this calculation 
is satisfactory and the maximum error, for N > 3, can be rated not to exceed 5%. Results 
and conclusions in complete agreement with ours are reached in the somewhat different 



analysis of the HT series presented in Ref. p3 . 



V. CONCLUSIONS 



We have presented our estimates of the low temperature quantities bo(N), c$(N) etc. de- 
fined by @, fllQ|) and (|TT|), obtained by a procedure which can essentially be seen as a simple 
improvement of the "matching method" introduced long ago in Ref. |35[ and since used sev- 



eral times with more or less unconvincing results, due either to inadequate implementation 
and/or to incorrect supplementary assumptions. The initial paper |)5| is an example of the 
former defect: the low temperature behavior was inadequately accounted for by one loop 
perturbation expansion and, on the HT side, too short series were used resulting into an 
unreliable matching. On the other hand Ref. is an example of both shortcomings since 
the use of HT series (at that time extending to ten terms only) was supplemented with the 
(now appearing obviously incorrect) conjecture that x an d A*2 have all positive HT coeffi- 
cients. Indeed even if we made the weaker assumption that there are at most finitely many 
negative expansion coefficients this would clearly imply that the nearest singularity of (for 
example) x is located on the real positive (3 axis. If also asymptotic freedom holds, then x 
should be a regular analytic function in the whole finite complex f3 plane, contrary to the 
numerical evidence presented in the previous section. 

We have tried to avoid the shortcomings of the previous approaches by the simplest pos- 
sible treatment of sufficiently long HT series and by excluding unwarranted supplementary 
assumptions. 

A brief review of some earlier references which are closely related to our analysis already 
appears in ||. Here we shall mention only some later studies and address the reader to Ref. 
for a long (but surely still incomplete) list of the abundant prior literature. 
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It is worthwhile to recall that recently, in the N = 3 case, a new method for extrapo- 
lating finite volume MC data to infinite volume [51] has been used to test the onset of the 
asymptotic behavior (P) by obtaining the second moment correlation length £ up to values 
as large as 10 5 lattice units and agreement with (H) has been found within 4% at (3 = 1. 
where £ ~ 10 5 . New MC data are now available ITD] also for the susceptibility which yield 



c (3) = 0.0146(10), cf c (4) = 0.0383(10) and c* fC (8) = 0.103(2) in very good agreement 
with our estimates c x (4) = 0.0344(3) and c x (8) = 0.1037(4) as well as with the conjectured 
exact results c x xact (A) = 0.0385 and c x xact (8) = 0.1027. Analogous results for the correla- 
tion length in the N = 3 case had been presented also in a computation extending to 
£ ~ 15000. In an earlier high precision multigrid MC study devoted to the N = 4 model 
on lattices of size up to 256 2 []5(J > the asymptotic behavior of £ and x na d been found to be 
perfectly compatible with (||) and ([H]) and the quantities <%(4) and c x (4) had been estimated 
to be c ? (4) = 0.0342(20) and c x (4) = 0.0329(16). Moreover in that study the possibility of 
an ordinary critical point (3 C < 1.25 was excluded, and it was stressed that the data could 
be compatible with a value f3 c > 1.25 only assuming implausibly large values for the critical 
exponents. ( A far away power singularity with a large exponent is likely to be merely an 
effective representation of an exponential behavior.) Also the MC single cluster simulation 
of Ref. [B7| for the N = 4 and N = 8 models gave good support to the asymptotic formulas 
CPTOp and (|9|) and produced estimates for completely consistent with (p"5|). 

Finally, on the side of the analytic approaches, we should mention the study of the scaling 
behavior in Refs. f68|| , whose results include a computation of the leading term of the 1/N 



expansion of c^, in complete agreement with the exact result (0), and of the same expansion 
for c x . 

In conclusion, we can summarize our main results as follows: 

a) By this and previous work || we have shown that our general N HT series are a useful 
tool also for obtaining high precision estimates of critical parameters in somewhat uncon- 
ventional contexts, giving further support to qualitative and quantitative results obtained 
by entirely different approximation methods. 

b) In the —2 < N < 2 vector models case we have confirmed, with high accuracy, 
the explicit formulas obtained by (semirigorous) model solving, for the critical exponents 
j(N), v(N) and A(iV). We have also computed the " dimensionless renormalized four point 
coupling constant" g r {N) for N = 0, 1,2 in complete agreement with other estimates, but 
with higher accuracy. 

c) For the N > 3 vector models, we can somehow extrapolate the HT series to the border 
of (or beyond) their region of convergence reliably enough to reconstruct the quantitative 
features of low temperature behavior and we can obtain a set of (hardly accidental) con- 
sistency checks with the predictions of the perturbative RG, of exact solutions and of 1/N 
expansions with an accuracy practically uniform with respect to N. As shown by Table 1, 
our estimates of the parameter c%(N) agree well with the exact calculation by the Bethe 
Ansatz (under the assumption (|i~5"D). On the other hand our estimates for c x (N), c r (N) and 
g r {N) are completely consistent with their 1/N expansions. 

Of course we must say that, strictly speaking, purely numerical computations cannot 
validate the RG predictions: only complete proofs can settle the question, but they are still 
to come. Therefore, in principle, further discussion of this subject may still be considered 
healthy and welcome as long as it may stimulate either to design a rigorous justification 
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of the generally accepted RG picture or to produce viable mechanisms for evading the 
expected asymptotic freedom regime while respecting the now established heuristic evidence. 
We believe, however, that the continuing accumulation of unambiguous, consistent and 
increasingly accurate numerical support for the RG predictions from a variety of independent 
approaches leaves little if any space for alternative pictures. 
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APPENDIX A: THE SUSCEPTIBILITY 

The HT coefficients of the susceptibility x{N,fl) = 1 + HT=\ o- r {N)l3 r are 
ai (N) = A/N 
a 2 (N) = 12/N 2 

a 3 (N) = (72 + 32iV)/(iV 3 (2 + N)) 

a 4 (N) = (200 + 76iV)/(iV 4 (2 + N)) 

a 5 (N) = 8(284 + 147iV + 20N 2 )/(N 5 (2 + N)(A + N)) 

a 6 (N) = 16(780 + 719iV + 201iV 2 + 19iV 3 )/(iV 6 (2 + iV) 2 (4 + N)) 

For the coefficients which follow it is typographically more convenient to set a r (N) = 
P r (N) I 'Q r (N) and to tabulate separately the numerator polynomial P r (N) and the denom- 
inator polynomial Q r (N). 

P 7 (N) = 16(26064 + 38076iV + 20742iV 2 + 5280iV 3 + 655iV 4 + 32iV 5 ) 

Q 7 (N) = N 7 {2 + N) 3 (A + N)(Q + N) 

P 8 (N) = 4(283968 + 383568iV + 186912iV 2 + 410007V 3 + 4392/V 4 + 187/V 5 ) 
Q 8 (N) = N 8 (2 + iV) 3 (4 + N)(6 + N) 

P 9 (N) = 8(3123456 + 4186336iV + 2087128A r2 + 492220iV 3 + 62386A^ 4 + 4161iV 5 + 116A^ 6 ) 
Qg(N) = 7V 9 (2 + N) 3 (A + N)(6 + N)(8 + N) 

P 10 (N) = 16(33868800 + 66758016/V + 53214272/V 2 + 22126648/V 3 + 5211372iV 4 + 
719330/V 5 + 58789/V 6 + 2684/V 7 + 55iV 8 ) 

Qio(N) = N 10 (2 + iV) 4 (4 + iV) 2 (6 + N)(8 + N) 

P n (N) = 32(3695370240 + 9913385984iV + 114372892 16 iV 2 + 7427564992iV 3 + 
2989987696iV 4 + 776848144iV 5 + 132130072iV 6 + 14693596iV 7 + 105291 IN 8 + 
46923iV 9 + 1225iV 10 + 16iy n ) 

Q n (N) = N u (2 + Nf(A + N) 3 (Q + N)(8 + N)(10 + N) 

P 12 (N) = 16(4990955520 + 11511967232iV + 10992991488iV 2 + 56098883527V 3 + 
1649559472iV 4 + 281912408/V 5 + 27080244/V 6 + 1334568/V 7 + 22368/V 8 - 199iV 9 + 5iV 10 ) 
Q 12 (N) = N 12 {2 + iV) 5 (4 + N) 2 (Q + N)(8 + N)(10 + N) 

P l3 (N) = 64(162478080000 + 406158981120iV + 431982472192iV 2 + 254291324928iV 3 + 
90288340864iV 4 + 19721001832iV 5 + 2561904944iV 6 + 170376718iV 7 + 1211742iV 8 
- 616479iV 9 - 37625iV 10 - 635iV n + 4iV 12 ) 

Q n (N) = N 13 (2 + iV) 5 (4 + N) 3 (6 + N)(8 + N)(10 + N)(12 + N) 
P U (N) = 16(21007560867840 + 63770201063424iV + 84400316350464iV 2 + 
63787725946880iV 3 + 30245054013440iV 4 + 9275137432448iV 5 + 
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1810742519232iV 6 + 2042842900 16 iV 7 + 7651128688iV 8 - 1135009056iV 9 - 177417600iV 10 

- 9861216iV 11 - 161554iV 12 + 5277iV 13 + 172iV 14 ) 

Q U (N) = iV 14 (2 + iV) 6 (4 + iV) 3 (6 + iV) 2 (8 + 7V)(10 + N)(12 + N) 
P 15 (N) = 16(9537349044142080 + 34641374485413888iV + 56169030631292928iV 2 + 
53537028436525056iV 3 + 33216830381735936iV 4 + 14006542675035136iV 5 + 
4047829991 104000iV 6 + 777531907925504iV 7 + 87227510881024iV 8 + 1944102682560iV 9 

- 1061882170400iV 10 - 183809104832iV 11 - 14563826832iV 12 - 515944376iV 13 + 
5830192JV 14 + 1259012iV 15 + 436477V 16 + 512/V 17 ) 

Q 15 (N) = 7V 15 (2 + 7V) 7 (4 + 7V) 3 (6 + iV) 3 (8 + N)(10 + N)(12 + N)(U + N) 
P 16 (N) = 4(410123375487221760 + 1537129944780374016iV+2583983411690471424iV 2 + 
2566975595695570944iV 3 + 1669084283351334912iV 4 + 7411 1401471 1103488iV 5 + 
225948162044579840iV 6 + 45385102417264640iV 7 + 4996850176026624iV 8 

- 64804204496896iV 9 - 122658733213440iV 10 - 20909429640960A^ n - 1752208241536iV 12 - 
56642417728iV 13 + 3062606512iV 14 + 412508368iV 15 + 18713696iV 16 + 

395328iV 17 + 3083iV 18 ) 

Q m (N) = N w (2 + iV) 7 (4 + iV) 4 (6 + iV) 3 (8 + N)(10 + N)(12 + N)(U + N) 
P 17 (N) = 8(35361815028050165760+ 138539666887258669056iV + 
245436909326998437888iV 2 + 259375081 142913859584iV 3 + 181396134616565809152iV 4 + 
87793764370648399872iV 5 + 29673202500166647808iV 6 + 6770762601709142016iV 7 + 
893508862130341888iV 8 + 5229394076767232iV 9 - 24934992828139008iV 10 

- 5547918408527104iV 11 - 625740097598720iV 12 - 33521607263744iV 13 + 
738107699392iV 14 + 272358030048iV 15 + 21867513640iV 16 + 937447020iV 17 + 
22261658iV 18 + 250495iV 19 + 692iV 20 ) 

Q 17 (N) = N 17 (2 + iV) 7 (4 + iV) 5 (6 + iV) 3 (8 + iV)(10 + N)(12 + iV)(14 + N)(16 + N) 
P 18 (N) = -16(-758936838540424642560- 3343934774878956158976iV- 
6720650800795024883712iV 2 - 8141819950133007089664iV 3 - 6611686534391523180544iV 4 
-3777832155850593533952iV 5 -1543669324445646061568iV 6 -443919373830353158144iV 7 - 
82653304109539049472iV 8 - 6333077582288386048iV 9 + 1419388196952978432iV 10 + 
578178922758906368iV 11 + 99290612095487744iV 12 + 9300735775467264iV 13 + 
256677768200576iV 14 - 53852331942080iV 15 - 8516631212960iV 16 - 629479458104iV 17 

- 26896421724iV 18 - 607483694iV 19 - 2912825iV 20 + 1540807V 21 + 2313iV 22 ) 

Q 18 (N) = N 18 (2 + iV) 8 (4 + N) 5 (6 + iV) 3 (8 + iV) 2 (10 + N)(12 + N)(U + N)(1Q + N) 
P 19 (N) = -32(-293792962132985669222400 - 1452203904509587992084480iV - 
3305764874023051160715264iV 2 - 45874982722 16547279765504iV 3 

- 4326547244550747303444480iV 4 - 2922303204727243671601152iV 5 

- 1446836388063996470624256iV 6 - 524786164553898279829504iV 7 

- 134502578329442459254784iV 8 - 21104878727348421885952iV 9 

-411991601001072488448iV 10 +794587338452494176256iV 11 +242141294836583751680iV 12 + 
39373307992978213888iV 13 + 3583854665917282560iV 14 + 51879072941552128iV 15 

- 36069375006840576iV 16 - 5868286096676352iV 17 - 508264525824336iV 18 

- 27200961065872iV 19 - 811699909040iV 20 - 3015005636iV 21 + 793163459iV 22 + 
32254806iV 23 + 562185iV 24 + 3824iV 25 ) 

Q 19 (N) = A^ 19 (2 + iV) 9 (4 + A^) 5 (6 + A^) 3 (8 + A^) 3 (10 + iV)(12 + A^)(14 + A^)(16 + A^)(18 + Ar) 
P 20 (N) = 25184031413058833177640960+ 120991848351738482367922176iV 
+ 2667357585644628251592622087V 2 + 356790558744797070187560960iV 3 + 
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322230995604339689974136832iV 4 

+ 206410863077042429645291520iV 5 + 95396352174319203631759360iV 6 + 
31338576528595789665009664iV 7 

+ 6741729364322236678275072iV 8 + 606442638553174662709248iV 9 
158967889827748034248704iV 10 

78160816104265974611968iV 11 - 16419585036974248984576iV 12 

1908913019215816540160iV 13 

- 62769211853172834304iV 14 + 19736004882625224704iV 15 + 4130688305419677696iV 16 
+ 421867767284303872iV 17 + 25241612021992960iV 18 + 693193236915968iV 19 

- 18814089206912iV 20 - 2684704080320iV 21 - 120121949760iV 22 - 2872757568iV 23 

- 35919232iV 24 - 178096iV 25 

Q 20 (N) = iV 20 (2 + A^) 9 (4 + Ar) 5 (6 + iV) 3 (8 + A^) 3 (10 + A^)(12 + A^)(14 + A^)(16 + A^)(18 + iV) 
P 21 (N) = 1351534773860942603511398400 + 6365460757030282723495772160iV 
+ 13733120620155454487896522752iV 2 + 17929816694573858749176348672iV 3 
+ 15741441712440400461107822592iV 4 + 9736866800557416285986619392iV 5 
+ 4291917210346656516848746496iV 6 + 1307842644179764469724872704iV 7 
+ 238056916142751738992001024iV 8 + 3543432139088928241090560iV 9 

- 12741 108714260109576372224iV 10 - 4295859766813135292465152iV 11 

- 755125851031606609051648iV 12 - 65023497347980126945280iV 13 
+ 2821567036764628910080iV 14 + 1727752072630205923328iV 15 

+ 264572837767576051712iV 16 + 22840865368902557696iV 17 
+ 1035703381014509568iV 18 - 6131507476388352iV 19 

- 4493143518510080iV 20 - 338088589058432iV 21 - 14045533700352iV 22 

- 355361402880iV 23 - 5200818400iV 24 - 35916176iV 25 - 47360iV 26 

Q 21 (N) = N 21 (2 + N) 9 (A + N) 5 (6 + N) 3 (8 + Nf(10 + N)(12 + N)(U + N)(lQ + N)(18 + 
N)(20 + N) 

In particular for N = we have (in terms of the variable (3 = (3/N): 
X(0, P) = 1 + 4/3 + 12/5 2 + 36/3 3 + 100/3 4 + 284/3 5 + 780^ 6 + 2172/3 7 + 5916/3 8 + 16268/3 9 + 
44100/3 10 + 

120292/3 n + 324932/3 12 + 881500/3 13 + 2374444/3 14 +6416596/3 15 + 17245332/3 16 +46466676/3 17 + 
124658732/3 18 + 

335116620/3 19 + 897697164/3 20 + 2408806028/3 21 . 
For N = 2 we have: 

X (2, (3) = 1 + 2(3 + 3(3 2 + 17/4/3 3 + 11/2/? 4 + 
329/48/3 5 + 529/64/3 6 + 14933/1536/5 7 + 5737/512/3 8 + 

389393/30720/? 9 + 2608499/184320/5 10 + 3834323/245760/5 11 + 1254799/73728/3 12 
+ 84375807/4587520/3 13 + 6511729891/330301440/3 14 + 66498259799/3170893824/5 15 + 
1054178743699/47563407360/3 16 + 39863505993331/1712282664960/5 17 + 

19830277603399/815372697600/3 18 + 

8656980509809027/342456532992000/3 19 + 2985467351081077/114152177664000/3 20 + 
811927408684296587/30136174903296000/3 21 
For N = 3 we have: 

X(3, f3) = 1 + 4/3(3 + A/3(3 2 + 56/45/3 3 + 428/405/? 4 + 
1448/1701/? 5 + 28048/42525/3 6 + 314288/637875/3 7 + 
685196/1913625/5 8 + 6845144/27064125/3 9 + 1159405664/6630710625/3 10 + 
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643017322016/5430552001875/3 11 + 915294455744/1 1636897146875/3 12 + 
12550612712128/244374840084375/? 13 + 120892276630256/3665622601265625/3 14 + 
3896992088570128/186946752664546875/3 15 + 
16982959056655084/1308627268651828125 / 3 16 + 
23336075876557256/2930091359102578125/3 17 + 
6224368647227625667744/1292302143475396569140625/3 18 + 
22688150130310720609472/7897401987905201255859375 / 3 19 + 
120072005214715268585744/71076617891146811302734375 ( 5 20 + 
43219200109596671558015312/44138579710402169818998046875/3 21 . 



APPENDIX B: THE SECOND CORRELATION MOMENT 

The HT coefficients of the second correlation moment fj, 2 (N,(3) = J2*?Li s r (N)(3 r are: 

Sl (N)=4/N 

s 2 (N) = 32 /N 2 

s 3 {n) = (328 + 16(W)/(iV 3 (2 + N)) 
s 4 (N) = (1408 + 64(W)/(iV 4 (2 + N)) 

s 5 (N) = (21728 + 14232N + 2208iV 2 )/(iV 5 (2 + iV)(4 + TV)) 
s 6 (N) = (156928 + 171072iV + 59840N 2 + 6784iV 3 )/(iV 6 (2 + iV) 2 (4 + N)) 
For the higher order coefficients s r (N) it is convenient to set s r (N) = P r (N)/Q r (N) and 
to tabulate separately the numerator P r {N) and the denominator Q r (N). 

P 7 (N) = 6487296 + 10904512iV + 7052384N 2 + 2192384N 3 + 328688iV 4 + 18944iV 5 
Q 7 (N) = N 8 (2 + N) 3 (A + N)(6 + N) 

P 8 (N) = 21596160 + 34468352iV + 21040128N 2 + 6162816iV 3 + 878208iV 4 + 48640iV 5 
Qs(N) = N 9 (2 + N) 3 (4 + N)(6 + N) 

P g (N) = 560007168 + 912207616iV + 585628864iV 2 + 190951904iV 3 + 33905168N 4 + 
3117448N 5 + 115616N 6 

Qg(N) = N 9 (2 + N) 3 (A + N)(6 + N)(8 + N) 

P 10 (N) = 14220853248 + 32437182464N + 31140450304iV 2 + 16449182208iV 3 + 
5251439360iV 4 + 1045045888N 5 + 127390272N 6 + 8712896iV 7 + 255616iV 8 
Q 10 (N) = N 10 (2 + iV) 4 (4 + iV) 2 (6 + N)(8 + N) 

P n (N) = 3549643407360 + 10748529770496iV + 14330317561856N 2 + 
11099228633088iV 3 + 5553033387520N 4 + 1888211571200N 5 + 
446700183296N 6 + 73788019072iV 7 + 8364424672iV 8 + 620464608N 9 + 
27094944iV 10 + 526848N 11 

Q n (N) = N n (2 + N) 5 (4 + N) 3 (Q + N)(8 + N)(1Q + N) 

P 12 (N) = 10920048721920 + 31897824264192iV + 
40878320844800N 2 + 30320450846720iV 3 + 14475616055296iV 4 + 4684012402688N 5 + 
1053225265152iV 6 + 165529253888N 7 + 17912071424N 8 + 1274573696iV 9 + 
53696896N 10 + 1013248iV n 

Qu(N) = N 12 (2 + iV) 5 (4 + N) 3 (Q + N)(8 + N)(10 + N) 

P 13 (N) = 398311403028480+ 1152646176964608N+ 1467863723606016N 2 + 
1086884435984384iV 3 + 521655212892160N 4 + 171508785347072iV 5 + 39822249767936iV 6 + 
6624256364416iV 7 + 788536224640 A^ 8 + 65778403648N 9 + 3655303744N 10 + 121407936N 11 + 
1818880N 12 
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Qi 3 (N) = N 13 (2 + iV) 5 (4 + N) 3 (Q + N)(8 + N)(10 + N)(12 + iV) 

P l4 (N) = 14381419069440000 + 49571261401006080iV + 76559245123780608iV 2 + 
70130377503539200iV 3 + 42548875902058496iV 4 + 18101 157266890752iV 5 + 
5581472628355072iV 6 + 1272436773582848iV 7 + 216726437780480iV 8 + 27616170103808iV 9 
+ 2605954492416iV 10 + 177259542528iV 11 + 8230177408iV 12 + 233423296iV 13 + 3045888iV 14 

Q 14 (jV) = iV 14 (2 + iV) 6 (4 + iV) 3 (6 + iV) 2 (8 + N)(10 + N)(12 + N) 

P 15 (N) = 7205314087322910720 + 29316772904369651712iV 
+ 54344822673767399424iV 2 + 60875749209938067456iV 3 + 46110757538166996992iV 4 + 
25056228677107990528iV 5 + 10119823516255248384iV 6 + 3107917444672675840iV 7 
+ 736807917027831808iV 8 + 1361445840304773 12iV 9 + 19700616020237824iV 10 + 
2231029399055360A^ n + 196196174024448iV 12 + 13160346699904iV 13 + 651384063232iV 14 + 
22420625856iV 15 + 478573168iV 16 + 4759552iV 17 

Q 15 (N) = N 15 (2 + N) 7 (A + N) 3 (Q + iV) 3 (8 + TV) (10 + iV)(12 + 7V)(14 + N) 

P ie (N) = 14215559923441336320 + 57338919520296763392iV + 
105267543605621293056iV 2 + 116619791028818280448iV 3 + 87189583554941026304iV 4 
+ 46640416019329581056iV 5 + 18481743700474920960iV 6 + 5546746731056758784iV 7 + 
1279666393060327424iV 8 + 229248311936622592iV 9 + 32092209212850176iV 10 + 
3517151947112448iV 11 + 300336273417216iV 12 + 19680328220160iV 13 + 958823550720iV 14 + 
32735074432iV 15 + 697558656iV 16 + 6957056iV 17 

g 16 (iV) = iV 16 (2 + iV) 7 (4 + iV) 4 (6 + iV) 2 (8 + N)(10 + N)(12 + JV)(14 + N) 

P 17 (N) = 16044762062435301457920 + 70470061100481487306752iV + 
142655241 108775004798976iV 2 + 176761774102693835440128iV 3 + 
150198190696943073624064iV 4 + 92960797041997201276928iV 5 + 
4346536 1502957674430464iV 6 + 15723355668874622926848iV 7 + 

4473127620779360665600iV 8 + 1012537718607308627968iV 9 + 183938871531499802624iV 10 + 
26986540252155746304iV 11 + 3211089966275868672iV 12 + 310291779078352896iV 13 + 
24273059658657280iV 14 + 1521460496772864iV 15 + 74832996108608iV 16 + 
2783265847904iV 17 + 73449080400iV 18 + 1221110008iV 19 + 9573792iV 20 

Q 17 (N) = N 17 (2 + iV) 7 (4 + iV) 5 (6 + iV) 3 (8 + N)(10 + N)(12 + N)(U + N)(16 + N) 
P 18 (N) = 750029539465981990010880 + 36711 18091050821760319488iV + 
8335725431444873454551040iV 2 + 11665373867250167725424640iV 3 + 
11276915336572927456837632iV 4 + 8001145549336838213206016iV 5 + 
4322734533275908527095808iV 6 + 1821577961516403890585600iV 7 + 
608650745059201940193280iV 8 + 163152472998781338124288iV 9 + 
35389047160612482121728iV 10 + 6253363024804357767168iV 11 + 
905249991372291489792iV 12 + 107885738708892270592iV 13 + 
10627036460386877440iV 14 + 866704728760623104iV 15 + 
58356430103347200iV 16 + 3208488174566912iV 17 + 140753928694016iV 18 + 
4731693506944iV 19 + 113947558464iV 20 + 1740125248iV 21 + 
12590720iV 22 

Q 18 (N) = N 18 (2 + iV) 8 (4 + iV) 5 (6 + iV) 3 (8 + iV) 2 (10 + N)(12 + JV)(14 + N)(16 + N) 
P 19 (N) = 627690146360050712617943040 + 3421601811563754334552326144iV + 
8712245471445327025300045824iV 2 + 13774632993059184247962599424iV 3 + 
15165775997343900382418436096iV 4 + 12361884192755007745538654208iV 5 + 
774431 1739948445340207677440iV 6 + 3821697127832819813536759808iV 7 + 
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1511074094374096761920684032TV 8 + 484537143774775297319108608TV 9 + 
127124521376337134571356160TV 10 + 27474245822576876554944512iV 11 + 
4918326246120771429597184TV 12 + 733003199414399719473152iV 13 + 
914282996271 76181227520TV 14 + 9600359324677991481344TV 15 + 

853907491123648274432TV 16 + 64659218756176879616TV 17 + 417187052022444288(W 18 + 
227737937559751168iV 19 + 10327555689201152TV 20 + 377014913696896iV 21 + 
10570062577056TV 22 + 212003127744TV 23 + 2690873184TV 24 + 16165376iV 25 

Q 19 (TV) = A^ 19 (2 + A^) 9 (4 + AT) 5 (6 + A^) 3 (8 + Ar) 3 (10 + A^)(12 + A^)(14 + A^)(16 + A^)(18 + A^) 
P 20 (7V) = 1815178094017325774709719040 + 96667604507856081268443709447V 
+ 24005039860199234508542705664iV 2 + 369399259770925560066169896967V 3 
+ 394915063820848500770466693127V 4 + 31171265041933192528668917760iV 5 
+ 18848811341049400534841163776iV 6 + 89441693578586454405823856647V 7 
+ 33853599986666151444993802247V 8 + 1033678218058491052840452096iV 9 
+ 2566517036975720475477934087V 10 + 521238144540532161039564807V 11 
+ 87017452575447588098539527V 12 + 12004983564719826766397447V 13 
+ 1378738461911714307112967V 14 + 133282419864415715000327V 15 
+ 11013708013947111997447V 16 + 791078082714905968647V 17 
+ 4985190283611369472iV 18 + 273323887975057408iV 19 
+ 12659942377698304iV 20 + 473471 116241408iV 21 + 135091633635847V 22 
+ 2729112569607V 23 + 34556842247V 24 + 205757447V 25 

Q 20 (N) = A^ 20 (2 + A^) 9 (4 + 7V) 5 (6 + 7V) 3 (8 + 7V) 3 (10 + 7V)(12 + 7V)(14 + 7V)(16 + 7V)(18 + 7V) 
P 2 i(N) = 104523299485865383076403609600 + 5486014563724569276060165734407V 
+ 13418599813193623120241995284487V 2 + 20324925293218392391657979904007V 3 
+ 21370019520816484485933438074887V 4 + 16572433317241480712315385937927V 5 
+ 9833396985154395402724955914247V 6 + 457150573406528768021251489792 * 7 
+ 1691739854246120103896217026567V 8 + 503678546409180584398511144967V 9 
+ 121507366326333039349957918727V 10 + 23864464633649983729873715207V 11 
+ 3830080098883358295958487047V 12 + 504556085413318840254791687V 13 
+ 55013729957017010858229767V 14 + 5050263468966785801748487V 15 

+ 402561333522263317135367V 16 + 29032280527863208755207V 17 + 

1948256629904685527047V 18 

+ 120024312967746739207V 19 + 6451766622729349127V 20 + 285190504156954887V 21 
+ 9837234706644487V 22 + 251744150118407V 23 + 4479171019847V 24 + 49479294567V 25 
+ 256263687V 26 

Q 21 (7V) = TV 21 (2 + 7V) 9 (4 + 7V) 5 (6 + 7V) 3 (8 + 7V) 3 (10 + TV) (12 + TV) (14 + TV) (16 + TV) (18 + 
TV) (20 + TV) 

In particular for TV = we have (in terms of the variable j3 = (3/N): 

fi 2 (0j) =4/3 + 32^ 2 + 164/3 3 + 
704/3 4 + 2716/3 5 + 9808/3 6 + 
33788/3 7 + 112480/3 8 + 364588/3 9 + 
1157296/3 10 + 3610884/3 11 + 11108448/3 12 + 
33765276/3 13 + 101594000/3 14 + 302977204/? 15 + 
896627936/3 16 + 2635423124/? 17 + 7699729296/3 18 + 
22374323436/3 19 + 64702914336/3 20 + 186289216332^ 21 . 

For TV = 2 we have: 
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/i 2 (2, (3) = 2(3 + 8(3 2 + 81/4/? 3 + 42(3 A + 3689/48/3 5 + 6193/48/3 6 + 
312149/1536/3 7 + 19499/64/3 8 + 13484753/30720/3 9 + 28201211/46080/5 10 + 
611969977/737280/3 11 + 101320493/92160/3 12 + 58900571047/41287680/3 13 + 
3336209179/1835008/5 14 + 1721567587879/754974720/3 15 + 
16763079262169/5945425920/3 16 + 58931 18865913171/1712282664960/3 17 + 
17775777329026559/4280706662400/3 18 + 1697692411053976387/342456532992000/3 19 + 
41816028466101527/7134511104000/3 20 + 206973837048951639371/30136174903296000/3 21 . 

For iV = 3 we have: 

^(3, (3) = 4/3(3 + 32/9/3 2 + 808/135/3 3 + 3328/405/3 4 + 
84296/8505/3 5 + 1391872/127575/3 6 + 21454864/1913625/3 7 + 
62634752/5740875/3 8 + 1923459304/189448875/5 9 + 
25854552704/2841733125/3 10 + 23813358832544/3016973334375/3 11 + 
180728998866176/27152760009375 / 3 12 + 
148615553292224/27152760009375 / 3 13 + 
161300027551 13536/3665622601265625/3 14 + 
647957301000434704/186946752664546875/3 15 + 
420771056234707712/157035272238219375/3 16 + 
15946537289290290889672/7832134202881191328125 / 3 17 + 
655875829145233998723328/430767381158465523046875/3 18 + 
26576776651881149990183488/23692205963715603767578125/3 19 + 
24824955666477074672626688/30461407667634347701171875 / 3 20 + 
318405686546338787648327888/544920737165458886654296875/3 21 . 



APPENDIX C: THE SECOND FIELD DERIVATIVE OF THE SUSCEPTIBILITY 

The HT coefficients of the second field derivative of the susceptibility X4(N,(3) = jj ( — 

2 + EZid r (N)[3 r ) are: 
di(N) = -32/ 'TV 

d 2 (N) = -8(64 + 35iV)/iV 2 (2 + TV) 
d 3 (N) = -256(12 + 7iV)/iV 3 (2 + TV) 

d 4 (N) = -8(15328 + 20856TV + 8952TV 2 + 1169TV 3 )/TV 4 (2 + TV) 2 (4 + TV) 
For the higher order coefficients d r (N) it is convenient to set d r (N) = P r (N)/Q r (N) and 
to tabulate separately the numerator P r (N) and the denominator Q r (N). 
P 5 (N) = -64(8568 + 11666TV + 5033TV 2 + 656iV 3 ) 
Q 5 (7V) = TV 5 (2 + 7V) 2 (4 + 7V) 

P 6 (N) = -32(848448 + 1708560iV + 1320504iV 2 + 483386iV 3 + 824927V 4 + 5229TV 5 ) 
Q 6 (TV) = TV 6 (2 + TV) 3 (4 + 7V)(6 + TV) 

P 7 (N) = -256(413760 + 8192487V + 6230787V 2 + 2245697V 3 + 377247V 4 + 23607V 5 ) 
Q 7 (7V) = 7V 7 (2 + 7V) 3 (4 + 7V)(6 + TV) 

P 8 (7V) = -8(3160154112 + 88808704007V + 105798502407V 2 + 69525082247V 3 + 
27453201927V 4 + 

6646220967V 5 + 961371847V 6 + 75911507V 7 + 2504377V 8 ) 
Q 8 (7V) = 7V 8 (2 + 7V) 4 (4 + 7V) 2 (6 + 7V)(8 + TV) 
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P 9 (jV) = -64(1423288320 + 3917961472iV + 4569267968iV 2 + 2938300752iV 3 + 
1135415360iV 4 + 269158808iV 5 + 38175104iV 6 + 2962181iV 7 + 96280iV 8 
Qg(N) = N 9 (2 + N) 4 (A + iV) 2 (6 + N)(8 + N) 

P 10 (N) = -32(797183508480 + 2819876773888iV + 4414993874944iV 2 + 
4031840569344iV 3 + 2382096324608iV 4 + 954513246240iV 5 + 
264332662544iV 6 + 50539833080iV 7 + 6533960024iV 8 + 543728966iV 9 
+ 26196458iV 10 + 553189iV n ) 

Q 10 (N) = N 10 (2 + iV) 5 (4 + N) 3 (Q + N)(8 + N)(10 + N) 
P n (N) = -256(340762460160 + 1180158980096iV+ 1807166870528iV 2 + 
1612510528512iV 3 + 9301 84344000 A^ 4 + 363788558192iV 5 + 
98342090128iV 6 + 18368180928iV 7 + 2323043092iV 8 + 189484599iV 9 + 
8969234iV 10 + 186536iy n ) 

Q n (N) = N u (2 + iV) 5 (4 + N) 3 (6 + N)(8 + N)(10 + N) 
P 12 (N) = -32(1315703150346240 + 54382178132951047V + 
10169455328329728/V 2 + 11390757326471168/V 3 + 85290383379660807V 4 + 
4512300585885440iV 5 + 1738422710398464/V 6 + 495320225789376iV 7 + 
104883376470528/V 8 + 16428534330480iV 9 + 1874379299088iV 10 + 
151083753396iV 11 + 8133230252iV 12 + 261683156/V 13 + 3795185iV 14 ) 
Q 12 (N) = 7V 12 (2 + iV) 6 (4 + iV) 3 (6 + 7V) 2 (8 + JV)(10 + iV)(12 + N) 
P 13 (N) = -512(270976473169920+ 1096725428969472/V + 
2005753217814528/V 2 + 2194539397009408/V 3 + 1603311500740608/V 4 + 
826889848699616/V 5 + 310370441728160iV 6 + 86142909672256/V 7 + 
17775181813000 iV 8 + 2715815783098/V 9 + 302703436000/V 10 + 
23883276726/V 11 + 1261330187/V 12 + 39906694/V 13 + 570432/V 14 ) 
Q 13 (N) = 7V 13 (2 + 7V) 6 (4 + iV) 3 (6 + N) 2 (8 + JV)(10 + N)(12 + N) 
P U (N) = -64(4729747131240284160 + 23387807447028596736/V + 
53391491 102489444352/V 2 + 74720767414227566592/V 3 + 
71798269494072377344/V 4 + 50281426991956893696/V 5 + 
26590502100001992704/V 6 + 108569394187044234247V 7 + 
3470103039376496128/V 8 + 874999045177634816/V 9 + 
1745557477264021 12/V 10 + 274995238990407047V 11 + 
3397617398241328/V 12 + 324891 1605006087V 13 
+ 23521314618812/V 14 + 1244052846554/V 15 + 45255618107iV 16 + 
1009996430/V 17 + 10401518iV 18 ) 

Q 14 (N) = iV 14 (2 + iV) 7 (4 + 7V) 4 (6 + iV) 3 (8 + JV)(10 + N)(12 + 7V)(14 + N) 
P 15 (N) = -512(1895033761431552000 + 9190498448855531520/V + 
20555882862533148672/V 2 + 28154602101819899904/V 3 + 
26448192225888174080iV 4 + 18088913492373442560/V 5 + 
9333532357592239 104iV 6 + 3715399595069458688iV 7 + 
1157147824100880000/V 8 + 284257647952320768/V 9 + 
55255671532695936/V 10 + 8487432924662416/V 11 + 
1023523926260184iV 12 + 95669381062712/V 13 + 6782408485072/V 14 + 
351973151712iV 15 + 12588857511iV 16 + 276796330/V 17 + 2813856iV 18 ) 

Q 15 (N) = iV 15 (2 + iV) 7 (4 + 7V) 4 (6 + Nf(8 + JV)(10 + iV)(12 + N) * (14 + TV) 
P 16 (N) = -8(393553748694956601507840 + 
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2238797016999938557476864iV + 5952616821086089377742848iV 2 + 
9835930646135145322512384iV 3 + 11329546912664804682366976iV 4 + 
9673321212788105516417024iV 5 + 6356189412895881926017024iV 6 + 
3294382649182089746055168iV 7 + 1369622314378844679962624iV 8 + 
462055197809790960812032iV 9 + 127470718364764484878336iV 10 + 
28892571739468927709184iV 11 + 5391343629949283622912iV 12 + 
827734092308104586240iV 13 + 104208728930679520256iV 14 + 
10686907762867520768iV 15 + 883185743012455424iV 16 + 
57860284710989984iV 17 + 2931823000694800iV 18 + 
110628321153680iV 19 + 2921354278696iV 20 + 48084058522iV 21 + 
370591 789 iV 22 ) 

Q 16 (N) = iV 16 (2 + iV) 8 (4 + iV) 5 (6 + iV) 3 (8 + iV) 2 (10 + N)(12 + JV)(14 + iV)(16 + N) 

P 17 (N) = -64(154328333132805848432640 + 
862171270522562364309504iV+ 2249151812525864836399104iV 2 + 
3642762103032724416626688iV 3 + 4108543937641919350308864iV 4 + 
3431288500060965921882112iV 5 + 2203091748524324305240064iV 6 + 
1114625140117096539160576iV 7 + 451935655809348535320576iV 8 + 
148577932801385167507456iV 9 + 399217396742292568391687V 10 + 
8810453666427880972288iV 11 + 1600832358513749923840iV 12 + 
239426089299260737792iV 13 + 29389521370216614400iV 14 + 
2942373642638070784iV 15 + 237767376253506048iV 16 + 
15259439999351088iV 17 + 758951651849152iV 18 + 28166343349968iV 19 + 
732948825144iV 20 + 11909447319iV 21 + 90757800iV 22 ) 

Q 17 (N) = N 17 {2 + iV) 8 (4 + iV) 5 (6 + iV) 3 (8 + iV) 2 (10 + N)(12 + N)(U + N)(16 + N) 

If we compute Nx 4 (N, (3) for N = O^we get: 

-2 - 32/3 - 256/3 2 - 1536/3 3 - 7664/3 4 - 34272/3 5 

- 141408/3 6 - 551680/3 7 - 2057392/3 8 - 7412960/3 9 - 25949984/3 10 

- 88740224/3 11 - 297422768/3 12 - 980094304/3 13 

- 3182108192/3 14 - 10199619200/3 15 

- 32321471824/3 16 - 101396444832/3 17 . 

For iV = 1 we have: 

X a(1,P) = -2 - 32(3 - 2Q4(3 2 - 4864/3/3 3 - 8232/3 4 - 553024/15/3 5 

- 2259616/15/3 6 - 180969728/315/3 7 - 217858792/105/3 8 - 20330135104/2835/3 9 

- 5377792736/225/3 10 - 12048694416128/155925/5 11 - 3450381618464/14175/3 12 

- 4559524221383168/6081075/3 13 - 32137492094329792/14189175/3 14 

- 4294238083842489856/638512875/3 15 - 66447301472480024/3378375/3 16 

- 615947855084824982464/10854718875/3 17 

For N = 2 we have: 

X4 (2, (3) = -1 - 8/3 - 67/2/5 2 - 104/3 3 - 12775 /48/3 4 - 1790/3/3 5 

- 931367/768/3 6 - 109691/48/3 7 - 93380347/23040/3 8 - 157557481/23040/3 9 

- 8158367639/737280/3 10 - 1061565359/61440/3 11 

- 1296061504531/49545216/3 12 - 477508721605/12386304/3 13 

- 439777014509471/7927234560/3 14 - 245861555567/3145728/3 15 

- 462463818305826161/4280706662400/3 16 - 314178739246240667/2140353331200/3 17 . 

For N = 3 we have: 
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X4 (3,/3) = -2/3 - 32/9/3 - 1352/135/3 2 - 2816/135/3 3 - 1520216/42525/3 4 

- 12992/243/? 5 - 3070624/42525/3 6 - 516883712/5740875/3 7 

- 697726412216/6630710625/3 8 - 331122359872/2841733125/3 9 

- 478066539947936/3878965715625/3 10 - 185574375218432/1481059636875/3 11 

- 150342773008769632/1221874200421875/3 12 - 429508071453349376/3665622601265625/3 

- 710293648879287815872/6543136343259140625/3 14 

- 1925804659821618529792/19629409029777421875/3 15 

- 2872493310184686424756616/33135952396805040234375/3 16 

- 32321239221821813512332352/430767381158465523046875/3 17 . 
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FIGURES 



Fig. 1. Numerical estimates of the critical exponent j{N) of the susceptibility, of the 
exponent v(N) of the correlation length and of the exponent 2A(N) from the second field 
derivative of the susceptibility, as computed for — 2 < iV < 2 by the method described in 
Section 3. Our results are represented by (the centers of) circles for the exponent v(N), 
squares for the exponent 'y(N) and triangles for twice the gap exponent A(iV) and they are 
compared with the corresponding exact formulas ([]) (dashed line) , (0) (continuous line) and 
(0) (dot-dashed line) respectively. Whenever no error bars appear, they are smaller than 
the data point. 



Fig. 2. Our estimates of the dimensionless renormalized coupling constant g r {N) for various 
values of < N < oo compared to some results from a recent MonteCarlo(MC) cluster 
computation for N = 2 and 3, to a field theoretic estimate in the case N = 3 and 
to other HT estimates [|4T1 , |3~2"f . For comparison, we have also plotted the large N asymptotic 
behavior (|2"7|). 

Fig. 3. The quantity A r (N) = (N/2) r a r (N) versus the order r for various fixed values 
of N. This normalization of the expansion coefficients a r (N) of the susceptibility has been 
chosen in order to make the plots for different values of N more easily comparable. We have 
also interpolated the data points by smooth continuous curves only to guide the eye. 



Fig. 4. The real part (circles) and the imaginary part (triangles) of the nearest singularity 
(3 n of the susceptibilty x(iV, (3) in the complex /3 = (3/N plane plotted as functions of 
x = 1 — 1/N in the range 1 < iV < oo. For x ^ 0.52, the singularity f3 n is still a real critical 
point, but for larger x it splits into a pair of complex conjugate singularities which move into 
the complex plane and, as iV — > oo, tend to the limiting points f3± ~ 0.32162(1 ± i). Here 
we have plotted only the trajectory of the singularity in the first quadrant of the complex {3 
plane. 



Fig. 5. The quantity B X (N, /3) defined by eq.flTTD versus (3 for N = 7 (lower set of curves) 
and N = 4 (upper set) showing the existence of a stationary point f3 s at which we estimate 
bo(N). We have plotted PA's of B X (N,{3) which use at least 15 HT series coefficients and 
with a difference between the degrees of numerator and denominator not larger than 4. 



Fig. 6. Numerical estimates of b (N) obtained starting from the quantities B X (N, (3 S ) (tri- 
angles), B^(N, $ S ) (circles) and B^N. ^s) (squares) are plotted versus N. Only for graphical 
convenience the estimates have been computed for three different sets of noninteger values 
of N. They are compared to the expected value (eq. ^) represented by the continuous line. 
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Fig. 7. Our numerical estimates of c%(N) (triangles), c x (N) (circles) and Cr(N) (squares) 
are compared to the exact fllBl) or conjectured |]6l| formulas (continuous lines), and to their 
1/N expansions (dashed lines). Please notice that for graphical convenience we have shifted 
upwards by 0.1 the data for c x (N) and have scaled down by a factor .25 the data for c r (N). 
The actual numerical values of these quantities are listed in Table 1. 
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TABLES 



TABLE I. Our central estimates of the universal constants c^(N), c x (N), c r (N) and of g r (N) 
for various values of N. We have indicated only the 'statistical' uncertainty. Beside our estimates 
of c^(N), c x (N), c r (N), we have reported in square parentheses the predicted value from eq. fll5|) 
or from the conjectured formula of Ref. |]6l| . 



N 




<bc( N ) 


C r (N) 


9r(N) 











10.53(2) 


1 








14.693(4) 


2 








18.3(2) 


3.5 


0.028(8) [0.0273] 


0.021(2) [0.0266] 


2.379(9) [2.45] 


20.4(1) 


4 


0.039(1) [0.0416] 


0.034(1) [0.0385] 


1.964(9) [2.007] 


20.9(1) 


5 


0.065(1) [0.0652] 


0.059(3) [0.0600] 


1.606(9) [1.624] 


21.6(2) 


6 


0.084(2) [0.0826] 


0.077(4) [0.0776] 


1.443(8) [1.452] 


21.9(3) 


8 


0.106(2) [0.1054] 


0.1035(8) [0.1027] 


1.290(9) [[1.291] 


22.5(4) 


10 


0.121(2) [0.1195] 


0.1212(5) [0.1194] 


1.213(5)[1.215] 


22.8(6) 


12 


0.130(2) [0.1290] 


0.134(2) [0.1312] 


1.169(5)[1.171] 


23.1(6) 


14 


0.137(2) [0.1358] 


0.143(1) [0.1399] 


1.141 (3) [1.141] 


23.3(6) 
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